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Report: ‘Ihe Carbon Cycle Revisited 

Bert Bolin and Inez Fung 

Introduction 

In developing an integrated perspective on the earth system, it has 
been essential to use simplified conceptual models representing the 
pool sizes and flux rates between major earth system components. The 
carbon cycle, because of the integral role it plays in the earth system, 
has served as a focal point for developing these models. From carbon- 
based models, new models have been developed, studies have been ini- 
tiated, and insights have been gained. To achieve greater understand- 
ing of the controls and dynamics of the earth system as represented by 
the conceptual model of the carbon cycle, it is beneficial to review our 
current state of knowledge regarding the carbon cycle to identify criti- 
cal gaps in our understanding of the earth system. 

Discussions during the Global Change Institute indicated a need 
to present, in some detail and as accurately as possible, our present 
knowledge about the carbon cycle, the uncertainties in this knowl- 
edge, and the reasons for these uncertainties. Figure 1 provides an 
overview of the carbon cycle. By its very nature, a diagram of this 
sort understates key assumptions and oversimplifies the Interlink- 
ing processes, so it may be misleading without rather precise anno- 
tations. Figure 1 is based on extensive research over a number of 
years and is not merely a rough qualitative overview, but all aspects 
of this research are not apparent from the condensed information in 
the figure. The extensive comments below offer a fuller appreciation 
of our present knowledge. Following these comments, we discuss 
basic issues of internal consistency within the carbon cycle, and 
end by summarizing the key unknowns. 
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Compartments 

Atmosphere 

Inventory: 748 ± 2 Gt C, based on the present (i.e., A.D. 1990) 
annual mean concentration of 353 ± 1 ppmv (1 ppmv COo = 2 12 
Gt C). 

Annual increase: Rate of 3.4 ± 0.4 Gt C corresponds to 1.6 ± 0.2 
ppmv/yr (cf. Watson et al., 1990; Tans et al., 1990). 

Oceans 

Total ocean volume : 1370 x 10 15 m 3 . The best estimate of the total 
amount of dissolved inorganic carbon (DIC) in the ocean is 38,000 ± 
700 Gt C, based on this volume and an average DIC concentration 
of 2.27 ± 0.04 mMol/kg. 

Surface layer: Usually defined as the top 100 m of warm surface 
water and about 200 m (sometimes more) of cold surface water. Vol- 
ume: ca. 50 x 10 15 m 3 ; the precise volume is somewhat arbitrary. 

Surface layer DIC: 1200 Gt; average concentration assumed to be 
2.02 + 0.05 mMol/kg (Takahashi et al., 1981); the uncertainty is 
partly dependent on the thickness of the layer. An annual increase 
of ocean surface carbon concentrations of about 0.05%/yr is taking 
place due to transfer of C0 2 from the atmosphere to the sea (assum- 
ing that quasi-equilibrium between the atmosphere and the sea is 
approximately being maintained and that the average Revelle factor 
is 10; see Rintoul report, this volume). Model calculations suggest 
that since preindustrial time an increase of this reservoir by 2-2.5%, 
i.e., 24-30 Gt C, has occurred (cf. Maier-Reimer and Hasselmann, 
1987; Bacastowand Maier-Reimer, 1990). 

Surface layer dissolved organic carbon (DOC): 50-250 Gt C; uncer- 
tainty due to few data available. Those available, if globally repre- 
sentative, show concentrations of about 0.2 mMol/kg, correspond- 
ing to a total amount of 120 Gt C. 

Surface layer biota: 3-5 Gt C, estimated indirectly by assuming 
that the amount present at any time is about 10% of the annual net 
primary production (NPP), i.e., the turnover time for biota in the 
surface layer is assumed to be about one month (de Vooys, 1979). 

Thermocline: On the average extending to about 1000 m in tropi- 
cal and midlatitude regions, absent in polar regions. Volume 300 x 
10 15 m 3 . 

Thermocline DIC: 8000 Gt C; average concentration assumed to be 
2-24 ± 0.06 mMol/kg (Takahashi et al., 1981); uncertainty is partly 
dependent on the definition of the size of the reservoir. This reser- 
voir is estimated to have increased by 75—125 Gt C since early last 
century due to transfer of excess carbon dioxide from the atmos- 
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phere directly or via the ocean surface layer (cf. Maier-Reimer and 
Hasselmann, 1987; Moore et al., 1989). 

Thermocline DOC: 300-500 Gt C, primarily in the upper half of 
the region; uncertainty due to very few data. The concentration 
seems to decline from the surface values given above to about 0.08 
mMol/kg in the lower part of the region. 

Deep sea (including the outcrop regions in the polar seas): Vol- 
ume 1020 x 10 15 m 3 . 

Deep sea DIC: 29,000 Gt C; average concentration assumed to be 
2.30 ± 0.03 mMol/kg (Takahashi et al., 1981); uncertainty is partly 
dependent on the definition of the size of the reservoir. 

Deep sea DOC: 700-1000 Gt C; average concentration assumed to 
be 0.06-0.09 mMol/kg; less relative uncertainty than for the surface 
layers because of a more homogeneous distribution as revealed by 
14c values of about 0.65, implying a radiocarbon age of 3000-4000 
years (cf. Williams, 1975). 

Terrestrial System 

Biota (nonwoody plus wood): 500-650 Gt C; assessment by direct 
inventory estimate. About 80% is judged to reside in wood (cf. Ajtay 
et al., 1979; Olson et al., 1983). 

Litter: 40-80 Gt C (cf. Ajtay et al., 1979; Olson et al., 1983); dead 
roots may not have been adequately accounted for. 

Soils: 1200-1500 Gt C (cf. Ajtay et al., 1979; Schlesinger, 1984); 
this estimate may to some extent include dead roots. It generally 
represents organic matter above about 50 cm depth. 

Peat: >150 Gt C (cf. Ajtay et al., 1979); this is probably an under- 
estimate, since peat often extends far below 50 cm, which usually is 
used for soil sample analysis. 

Fossil fuels: Estimates of fossil fuel resources are given in Table 1. 


Fluxes 

Atmosphere-Ocean 

This flux is 60-90 Gt C/yr. In steady state there is a mutual mol- 
ecular exchange maintained by the partial pressure of carbon diox- 
ide (p0O 2 ), preindustrially about 280 ppmv, in the atmosphere and 
the ocean. Given the 1973 pC0 2 of 328 ppmv, this yields an esti- 
mated average residence time for a C0 2 molecule in the atmosphere 
of 8-12 years and an estimated exchange rate of 15-21 mol 
C0 2 /yr/m 2 . These estimates (cf. Bolin et al., 1981) are based on (1) 
determination of rate of radon evasion from the sea (Peng et al., 
1979), (2) balance between 14 C radiocarbon invasion into the sea 
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Table 1: Fossil fuel reserves 


Coal 
(10 9 ton) 


Oil 

(10 9 m3) 


Gas 

(1012 m3) 


Used until 1988 

Proven reserves 
(1985) 

Probable reserves, 
low/mean /high* 


150 

700 

>5000 


100 

100 

55/135/275 


40 

95 

55/150/330 


1.03 tons of oil equivalent; 1 ton coal = 0.65 


*95/50/5% probability 
Energy content: 1000 m 3 gas = 
tons of oil equivalent. 

Carbon content: Coal and lignite, 69%; oil, 84%; gas, 540 g/m 3 
Relative C0 2 emission per unit of energy: gas, 1.0; oil, 1.5; coal 1 8 
Fr °988 Chilhng and WiCgand ' 1987: Whlte - 1987; cf - 3150 Goldemberg et al„ 


and radioactive decay in the sea, and (3) rate of invasion of 14 C pro- 
duced by nuclear bomb testing (Stuiver, 1980). 

During preindustrial times a net evasion of C0 2 from the oceans 
in tropical latitudes was approximately balanced by a net invasion 
into the oceans in polar and subpolar latitudes, which were associ- 
ated with poleward transfer in the atmosphere and equatorward 
transfer in the oceans. The net circulation is estimated to have been 
2-4 Gt C/yr (Pearman and Hyson, 1986; Keeling et al., 1989a, 
1989b). The derivation of this transfer rate is based on an evalua- 
tion of the transport capability of the atmosphere between tropical 
latitudes and polar regions and an assessment of differences 
between the average atmospheric concentrations of C0 2 in these 
regions of 1-1.5 ppmv, as evaluated from data taken around 1960 
(cf. Keeling et al., 1989a, 1989b). This pattern has been markedly 
disturbed due to anthropogenic C0 2 emissions, primarily in the 
Northern Hemisphere. Fossil fuel burning, deforestation, and chang- 
ing land use have created a mean excess of 5-10 ppmv of C0 2 in the 
atmosphere relative to the oceans. As a result, a net flux of C0 2 into 
the- oceans takes place at present, which is estimated to be 08-2.4 
Gt C/yr (cf. Bolin, 1986; Siegenthaler and Oeschger, 1987; Tans et 
al., 1990; Watson et al., 1990). 

Within the Oceans 

The following estimates are based both on direct measurements 
and on assessments with the aid of carbon cycle models (cf. 
Oeschger et al., 1975; Siegenthaler, 1983; Bolin, 1986; Maier- 
Reimer and Hasselmann, 1987; Sarmiento et al., 1992). 
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DIC in surface water used for marine biota production: 30-50 Gt 
C/yr; estimated from a large number of measurements of NPP. 

Marine biota transferred into dead organic matter, detritus, and 
DOC, and partly mineralized in the surface layers within a month to 
a year: 30-50 Gt C/yr. This flux is, however, less than the total NPP 
referred to above because a part is transferred into the thermocline 
region (see below). 

Detritus settling into the thermocline region: 2-6 Gt C/yr. Mineral- 
ization of detritus and DOC in the thermocline region: 2-6 Gt C/yr 
as estimated from the oxygen deficit (AOU, apparent oxygen utiliza- 
tion) in the layer. 

Transfer of detritus and DOC into the deep sea and mineralization 
into DIC: 1-3 Gt C/yr, estimated on the basis of radiocarbon age of 

DIC in the deep sea (see below). 

DIC (gross) exchange between the surface layer, thermocline 
region, and deep sea due to ocean circulation and turbulence. 25—40 
Gt C/yr, estimated with the aid of present knowledge of ocean circu- 
lation: of turbulent transfer, including penetrative convection in 
polar regions; and of the radiocarbon age of DIC in the deep sea. 
Due to a DIC gradient between the deep sea and the surface layer, 
this leads to: 

Net upward transfer of DIC to the thermocline region: 1-3 Gt 
C/yr, and to the surface layer: 2-6 Gt C/yr. Because of increased 
DIC in the surface layer during the 20th century caused by the 
transfer of excess carbon from the atmosphere to the sea, this 
upward transfer has decreased by 1-2.5 Gt C/yr, primarily from the 
thermocline region to the surface layer, which enhances the uptake 
by the oceans of excess atmospheric C0 2 (cf. above). Mineralization 
of detritus and/or DOC in the deep sea: 1-3 Gt C/yr, balancing the 
net upward transfer due to water motions (see above). 

The fluxes as given above constitute the biological pump, which 
maintains a higher pC0 2 (400-700 ppmv) in the thermocline region 
and the deep sea than would be expected otherwise. 


Atmosphere-Terrestrial Exchange 
(Preindustrial Estimates) 

The following are direct estimates of undisturbed preindustrial 
conditions as obtained from field data gathered in major biomes (cf. 
Ajtay et al., 1979; Olson et al., 1983; Bolin, 1986). 

• Gross primaiy production: 90-130 Gt C/yr 

. Respiration: 40-60 Gt C/yr, i.e„ NPP: 40-70 Gt/yr, of which 
about 25% is stored in wood 
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• Formation of dead organic matter, litter: 40-70 Gt C /yr 

• Litter decay within a few years: 35-65 Gt C /yr 

• Litter to organic carbon in soil: 2-6 Gt C/yr 

• Soil decomposition: 2-6 Gt C/yr 

• Peat formation as well as peat decomposition: <1 Gt C/yr. 


Erosion of the Continents and Flux to the Ocean 
(Preindustrial Estimates) 

The following estimates of preindustrial river loads have been cal- 
culated (cf. Kempe, 1979): 

• HC0 3 -C: 0.4 Gt C/yr 

• DOC: 0.1 Gt C/yr 

• Particulate organic carbon (POC): 0.1 Gt C/yr 

• Inorganic, suspended C: 0.2 Gt C/yr 

• Total: 0.8 Gt C/yr. 

It is likely that the flux of carbon to the sea has increased due to 
anthropogenic activities, this being particularly so for the flux of 
DOC and POC. 


Other Fluxes 

Carbon Jlux from the ocean to the sea Jloor: Sediments of carbon- 
ates: 0. 1-0.2 Gt C/yr; organic carbon, open sea: about 0.05 Gt 
C/yr; shallow seas, partially anthropogenic: 0.3-0. 5 Gt C/yr (cf 
Kempe, 1979). 

Deforestation and changing land use: Release of C0 2 into the 
atmosphere: 0.4-2. 6 Gt C/yr, accumulated release from about 1750 
to 1980: 80-140 Gt C. These direct estimates have been made by 
using information from national forest services, calibrated against 
field measurements from present activities as well as satellite obser- 
vations (Houghton et al„ 1987; Detweiler and Hall, 1988). 

Fossil fuel combustion: At present 5.6-6.0 Gt C/yr; accumulated 
emissions since the beginning of the industrial revolution (about 
1840) until 1985: 190-210 Gt C. At present, the respective contribu- 
tions from burning coal, oil, and gas are approximately 40%, 44%, 
and 16%. 


p wi in” 


Modeling the ‘Earth System 


158 


Constraints on the System 

Due to a number of constraints on the system, it is not permissi- 
ble to combine any numbers given for reservoir sizes and fluxes 
within the uncertainty ranges that have been assigned. These con- 
straints are listed below. 

. It should be noted that variations of atmospheric C0 2 . during the 
period from A.D. 1000 until anthropogenic Influences became sig- 
nificant, have been less than 3-5 ppmv for any period of a decade 
or two. Therefore, the mean net exchange between the atmos- 
phere, on the one hand, and the oceans and the terrestrial sys- 
tem on the other, must have been less than 1 Gt C/yr and prob- 
ably less than 0.5 Gt C/yr. This close balance is remarkable in 
the light of the sizes of the different reservoirs and the gross 
exchanges that take place between them. 

• This balance has been disturbed by emissions of C0 2 from the 
burning of fossil fuels, deforestation, and changing land use. The 
average fossil fuel input during 1980-1987 was 5. 1-5.5 Gt C/yr, 
while that from deforestation and changing land use was 0.4-2. 6 
Gt C/yr, i.e., the total anthropogenic emissions were 5.9-8. 1 Gt 
C/yr. The atmospheric increase during this same period was 
2. 8-3. 4 Gt C/yr, while the oceanic uptake is computed to have 
been 1. 6-3.0 Gt C/yr, i.e., the total sinks are 4.4-6.4 Gt C/yr. 
Although there is a slight overlap between the estimated ranges 
for the total sources and the total sinks, it seems likely that there 
is an additional unknown sink, the “missing sink." 

• In addition to showing an upward trend, the atmospheric C0 2 
data from the globally distributed network of monitoring sites also 
show that C0 2 concentrations are higher in the Northern Hemi- 
sphere than in the south, reflecting the dominance of the North- 
ern Hemisphere midlatitude fossil fuel sources. In 1980-87 the 
average concentration difference between 82°N (Alert) and the 
South Pole was 3.0 ppmv. Tans et al. (1990) attempted to simu- 
lated the north-south C0 2 gradient in a three-dimensional atmos- 
pheric tracer transport model. If we accept the geographic distrib- 
ution of the fossil fuel sources and an approximate distribution of 
land sources, a major sink is required in the Northern Hemi- 
sphere rather than in the Southern Hemisphere to match the 
atmospheric observations. 

• The difference in C0 2 concentrations between the hemispheres 
during the preindustrial era is unknown and is a source of uncer- 
tainty in the above inference about the hemispheric distribution 
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of the sinks. Keeling et al. (1989a, 1989b) have extrapolated the 
concentrations back in time and estimated that the C0 2 concen- 
tration at the South Pole may have been about 2 ppmv higher 
than in the Northern Hemisphere in the 1700s. Even when this 
estimate is included in the calculations, the requirement of a 
comparatively large contemporary northern sink remains. 

The data from the Northern Hemisphere further indicate that the 
magnitude of the sinks in the North Atlantic and the North Pacific 
was significantly less than what is required to balance the carbon 
budget and that a land sink for anthropogenic C0 2 must be 
assumed. The mechanism of this missing sink is unknown. 

A major uncertainty remains, however, in that the estimate of the 
uptake by the oceans critically depends on the determination of 
the global distribution of pC0 2 in ocean surface water. Observa- 
tions, particularly from the southern seas, are still few. 

Carbon cycles should be tested against data on total carbon, as 
well as against the distributions of 8 13 C and A 14 C. The transient 
changes of the latter due to anthropogenic activities also put 
some constraints on the behavior of the system as a whole. The 
decrease of A^C from the middle of last century until 1954 (the 
Suess effect) is thus consistent with the estimated past emissions 
of fossil carbon (free of 14 C), although the associated decrease of 
the concentrations in the surface water of the oceans has not 
been verified by measurements. 

M. Heimann and R. Keeling have carried out simulations using a 
three-dimensional ocean model, similar to those reported by Tans 
et al. (1990), including a simulation of the 13 C distribution in the 
1980s. They found that a large land sink in the Northern Hemi- 
sphere does not match the observed 8 13 C gradient. They confirm, 
however, that uptake by the Northern Hemisphere oceans would 
require a partial pressure difference between the atmosphere and 
the ocean far In excess of what is observed. A small land sink, 
approximately proportional to the net primary production, is in 
any case required to match available observations. 

It is further important to consider the history of the C0 2 balance 
since the 1700s and, in particular, to assess the integrated 
changes during this period. Houghton (1989) estimated both the 
total magnitude and the distribution with time of a missing sink 
for the period 1860-1985. Although such a residual is necessarily 
quite uncertain, it appears that there was no missing sink before 
about 1940. It was largest in the 1970s and has decreased since 
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then A careful analysis of such modeling results with gradually 
improved models might reveal the nature of this sink or the inad- 
equacy of other features of present carbon cycle models. 

Key Questions and Suggestions for Further Research 

The global carbon cycle remains a puzzle. The enigma is caused 
by the difficulty of making direct measurements and the paucity of 
such measurements. Increasingly complex models can be developed 
in order to estimate the relative importance of difference processes 
and interactions and to assess the sensitivity of atmospheric C0 2 
concentrations to assumed changes in these processes. Critical 
areas with remaining questions and important research needs are 
listed below. 

• What is the current distribution of pC0 2 in the oceans? What are 
the seasonal variations of oceanic pC0 2 ? What is the magnitude 
of the errors introduced by inadequate seasonal sampling in cur 
rent estimates of the annual mean pC0 2 , especially in the North 
Pacific and North Atlantic? What is the role of the marginal seas? 

• What are the rate-limiting processes that determine the magni- 
tude of the ocean uptake of anthropogenic C0 2 ? How can we 
improve our modeling of these processes? How can we improve 
the model calculations on time scales important for anthro- 
pogenic 0O 2 ? Currently, various tracers whose input histories dif- 
fer from that of fossil fuel C0 2 are used to calibrate/validate 
transport characteristics of various time scales in ocean models. 
The magnitude of the C0 2 uptake by the oceans in these models 
is very sensitive to the calibration procedure (e.g., Joos and 
Siegenthaler, 1989). 

• What was the preindustrial distribution of C0 2 in the atmosphere 
and in the oceans? Was the Southern Hemisphere concentration 
higher than in the north, as suggested by Keeling’s extrapolation 
of the Mauna Loa and South Pole records? If so, what processes 
were responsible for the elevated pC0 2 (relative to the atmos- 
phere) in the southern oceans? Simulation of preindustrial C 
distribution (Levin et al„ 1987) may help test hypotheses of air- 
sea C0 2 exchange. 

• Seasonal transport in atmospheric models of the global carbon 
cycle should be tested by the simulation of inert tracers. The vali- 
dation has only been done in the annual mean. 

• Establishment of a calibration standard for 13 C measurements 
and expansion of the network for monitoring atmospheric C 
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must receive high priority, so that latitudinal gradients and long- 
term trends can be interpreted with confidence. Likewise, carbon 
cycle models must include 13 C together with C0 2 to constrain the 
partitioning between land and ocean sources and sinks. 

• The interannual variations of C0 2 in the atmosphere provide 
clues to the carbon cycle’s responses to perturbations. What tran- 
sient shifts in the C0 2 balance between photosynthesis and 
decomposition have resulted from increasing C0 2 concentration 
and fluctuations in surface air temperature and precipitation? 
Similarly, the upper ocean temperature is known to have fluctu- 
ated in the past 30 years. How might these fluctuations in ther- 
modynamics and associated dynamics affect the oceanic carbon 
cycle and be manifested in the atmospheric C0 2 record? 

• How have the numerous carbon pools in the terrestrial biosphere 
been altered by deforestation, afforestation, and changing agricul- 
tural practices? How can we narrow the uncertainties in our esti- 
mates of the C0 2 sources that result from land use modification? 

The size of the different carbon pools and the dominant controls 
and turnover rates of each pool need to be determined. How 
would climate change and anthropogenic perturbations change 
substrate properties (such as pH and soil fertility) and alter pho- 
tosynthesis, allocation, decomposition, and other dynamical prop- 
erties of terrestrial ecosystems? What controls the geographic 
variations in these processes? Can the uncertainties of C flux into 
various pools be reduced to determine whether a Northern Hemi- 
sphere terrestrial sink for C exists? 
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